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ABSTRACT 

We cross-correlate the COB^ DMR 2-year sky maps with spatial templates from 
long- wavelength radio surveys and the far-infrared COBE DIRBE maps. We place an 
upper limit on the spectral index of synchrotron radiation /Sgynch < —2.9 between 408 
MHz and 31.5 GHz. We obtain a statistically significant cross-correlation with the 
DIRBE maps whose dependence on the DMR frequencies indicates a superposition of 
dust and free-free emission. The high-latitude dust emission {\b\ > 30°) is well fitted 
by a single dust component with temperature T = 181^7 K and emissivity e oc {p /pqY 
with = l.^tot The free-free emission is spatially correlated with the dust on 
angular scales larger than the 7° DMR beam, with rms variations 5.3 ± 1.8 /iK 
at 53 GHz and angular power spectrum P oc If this correlation persists to 

smaller angular scales, free-free emission should not be a significant contaminant 
to measurements of the cosmic microwave anisotropy at degree angular scales for 
frequencies above 20 GHz. 

Subject headings: radio continuum: interstellar - ISM: dust - cosmic microwave 
background 
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1. Introduction 



The microwave sky is dominated by the cosmic microwave background (CMB) 
and Galactic synchrotron, dust, and free- free emission. A number of authors 
have attempted to distinguish these components based on sky surveys where one 
component dominates the others and on the frequency spectrum inferred from 
microphysical processes (Fixsen, Cheng, & Wilkinson 1983; Lubin et al. 1985; Wright 
et al. 1991; Bennett et al. 1992; Bensadoun et al. 1993; Gutierrez de la Cruz et 
al. 1995). This approach has been most successful for synchrotron radiation: radio 
surveys and local measurements of the cosmic ray electron energy spectrum determine 
both the spectrum and morphology within broad limits (Banday & Wolfendale 
1991, Bennett et al. 1992). Free-free emission from electron-ion interactions has 
well-determined spectral behavior but is poorly mapped at high Galactic latitudes 
(Bennett et al. 1992, Reynolds 1992). Dust emission suffers the opposite problem: 
although it dominates infrared surveys, its spectral behavior depends on the shape, 
composition, and size distribution of the dust grains, which are poorly known (Desert, 
Boulanger, & Puget 1990). 

At millimeter wavelengths, the large angular scale anisotropy in the CMB 
is larger than fluctuations in the combined Galactic foregrounds, so that detailed 
Galactic modelling is not required for cosmological tests (Bennett et al. 1992, Smoot 
et al. 1992). At finer angular resolution or across a broader frequency band, this is not 
necessarily true (Banday & Wolfendale 1990, 1991; Franceschini et al. 1989). Precise 
maps of Galactic free-free and dust emission provide information on physical processes 
in the interstellar medium and aid planning of future CMB anisotropy measurements. 

Synchrotron emission results from the acceleration of cosmic-ray electrons in 
the Galactic magnetic field, and may be approximated as a power law in frequency, 
7"synch oc . Bennett et al. (1992) review several models of synchrotron emission based 
on radio surveys at 408 MHz (Haslam et al. 1981) and 1420 MHz (Reich & Reich 
1988). They consider three models for the spectral index: spatially invariant, spatially 
varying but frequency-independent, and a spatially varying index that steepens with 
frequency as determined by the local cosmic-ray energy spectrum. Comparison with 
a survey at 19 GHz (Boughn et al. 1992) suggests that the third model is somewhat 
better than the other two, although all have obvious shortcomings. Synchrotron 
indices used by various authors to scale the 408 MHz survey to millimeter wavelengths 
range from -2.7 to -3.1, corresponding to a factor of 7 uncertainty in amplitude at 53 
GHz. 

Bennett et al. (1992, 1994a) provide a full-sky map of the Galactic free- 
free emission derived from linear combinations of the COBE Differential Microwave 
Radiometers (DMR) sky maps, designed to isolate emission with Tg oc 
Although in principle this map presents an unbiased full-sky template for emission 
from the ionized gas in the interstellar medium, in practice the noise levels are such 
that only the quadrupolar component AT oc esc |6| is detected at |6|>15°. Other 
tracers of the warm ionized medium suffer from signal to noise ratio (e.g.. Nil emission. 
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Bennett et al. 1994b) or from undersampling and selection biases (e.g., Ha emission 

and pulsar dispersion measure, Reynolds 1992, 1984). 

Wright et al. (1991) provide estimates of far-infrared Galactic dust emission 
using the COBEYax Infrared Spectrophotometer (FIRAS). Assuming that the dust 
intensity along each line of sight may be characterized by optically thin emission with 
frequency-dependent emissivity 1^ = tq (^)'^ B,y(T) where B^, is the Planck function 
and z/q = 900 GHz, they find temperature T = 23.3 K and emissivity j3 = 1.65 
for the mean spectrum including the Galactic plane. A fit with (5 fixed at 2 shows 
excess emission at long wavelengths, characterized by a two-temperature model with 
Ti 20.4 K and T2 « 4.8 K. Reach et al. (1995a) extend this model with an improved 
data set; at high latitudes (|6|>30°) they detect a pervasive long wavelength excess 
characterized by either a flattened emissivity (3 = 1.2 ±0.3 or a combination of warm 
and cold components with Ti = 17.6 ± 0.5 K and T2 = 6.9 ± 0.8 K. The opacities of 
the warm and cold components are highly correlated; however, the FIRAS data do 
not distinguish between two distinct (although spatially correlated) dust populations 
or a single dust population with enhanced submillimeter emissivity. 

Attempts to measure Galactic emission at millimeter wavelengths are hindered 
by the small amplitude of the individual emission components and the presence both 
of other Galactic sources and of anisotropy in the cosmic microwave background. 
Point-by-point decomposition based on the different frequency dependences of the 
emission mechanisms are typically limited by signal to noise considerations (Brandt 
et al. 1994). The abihty to detect weak signals is significantly enhanced if the angular 
distribution is known a priori. In this paper, we cross-correlate the DMR full-sky 
maps with spatial templates from long-wavelength radio surveys and the far-infrared 
maps of the COBE Diffuse Infrared Background Experiment (DIRBE, Boggess et 
al. 1992). After accounting for instrument noise and chance alignments of the CMB 
anisotropy, we derive upper limits to the synchrotron spectral index steeper than 
used by many authors. We obtain a statistically significant cross-correlation with the 
DIRBE maps whose dependence on the DMR frequencies indicates a superposition 
of dust and free-free emission. The dust temperature and emissivity derived from 
angular variations in the DMR and DIRBE data agree with the values derived 
independently from the absolute FIRAS spectra, and place a lower hmit on enhanced 
dust emissivity at long wavelengths. The detection of free-free emission correlated 
with the DIRBE dust emission provides a spatial template and normalization for at 
least one component of the warm ionized medium at various angular scales. 
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2. Analysis 



We assume that the DMR microwave maps may be represented by a superposition 
of CMB anisotropy and Galactic signals whose angular distribution is traced by an 
external map, 

where AT^'^^ is the antenna temperaturei in a DMR map (in units mK), AX*^^' is 
the intensity of the Galactic template map (not necessarily in temperature units), and 
the coefficient a converts the units of the Galactic map to antenna temperature at 
the DMR frequency. Fitting for a on a pixel- by-pixel basis (e.g., a linear correlation 
coefficient) does not make full use of the angular information in the two maps; instead, 
we compare the full cross-correlation function 

c{e) = J2 ^x^"' ATf 

id 

(summed over all pixel pairs {i,j} separated by angle 9) to the autocorrelation 
function 

A{e) = AXf AXf ^1 

and find the coefficient a by minimizing 

= E(C - aAUM-')at,iC - aA), 

a,b 

where M is the covariance matrix of the cross-correlation function, and the indices 
a and b run over the angular bins of the correlation functions. Bennett et al. (1993) 
use the same technique to search for correlation between the DMR maps and various 
extragalactic sources. 

We use the 2-year DMR sky maps (Bennett et al. 1994a) at 9.5, 5.7, and 3.3 
mm wavelength (31.5, 53, and 90 GHz). For long wavelength (synchrotron) Galactic 
templates, we use either the 408 MHz survey (Haslam et al. 1981) or the 408 MHz 
survey scaled to the DMR frequencies using a spatially varying frequency-dependent 
spectral index as described in Bennett et al. (1992), denoted here as the "cosmic 
ray" model. For short wavelength (dust) Galactic templates, we use the DIRBE sky 
maps at 240, 140, and 100 /im wavelength from which a model of the zodiacal dust 
emission has been removed (Reach et al. 1995b). We convolve the Galactic template 
maps with the DMR beam pattern (Wright et al. 1994) and pixelize them in the same 
Galactic representation used for the DMR maps. 

The Galactic plane has considerable structure smaller than the 7° DMR beam 
size; we exclude the Galactic plane (|6|<20° or |6|<30°) from all analysis. The 
DMR maps are insensitive to a monopole and have a dipole dominated by the CMB. 

^ Antenna temperature is defined in terms of the power received per unit bandwidtli, 

2 

P = kTA^L' wliere k is Boltzmann's constant. It is related to tlic intensity Ii, by Ii, = 2kTA^- 
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Galactic microwave emission from synchrotron, free-free, and dust is dominated by 
a csc|6| quadrupole. Analysis which includes the quadrupole could result in a false 
positive cross-correlation between the DMR maps and a Galactic template due to the 
microwave emission from a different template with a brighter quadrupole at millimeter 
wavelengths. Furthermore, the GMB quadrupole has a component counter- aligned 
with a Galactic cosecant of similar magnitude (Bennett et al. 1992), potentially 
masking correlation between the DMR and Galactic maps. Accordingly, we remove 
the monopole, dipole, and quadrupole from the the high-latitude portion of both the 
DMR and Galactic template maps prior to computing the auto- and cross-correlation 
functions. As a consequence, our results depend only on the intensity variations in 
the maps and are insensitive to the zero level of either the DMR or Galactic template 
maps. 

We determine the covariance matrix M of the cross-correlation function and 
assess the statistical significance of the coefficient a. as Galactic features occasionally 
align with GMB anisotropy by using Monte Carlo simulations in which the DMR 
map is replaced by a scale-invariant Harrison-Zel'dovich spectrum of GMB anisotropy 
normalized to the quadrupole value Qrms-ps = 20 fiK determined from the 2-year 
DMR data (Gorski et al. 1994). To each GMB realization we add a reahzation of 
instrument noise defined by the level and pattern of noise in the DMR A and B 
channels, then combine the channels to form the (A+B)/2 sum map. We excise 
the Galactic plane and remove the monopole, dipole, and quadrupole from each 
simulated sum map prior to computing the cross-correlation. The mean coupling 
between angular bins defines the covariance matrix, M^b — j^Yl CaCb averaged over 
N simulations; throughout this work we use 72 angular bins each of width 2?6. The 
simulations thus automatically account for cosmic variance from random alignments 
of GMB anisotropy with features in the Galactic template maps as well as the effect 
of instrument noise and the aliasing of power between multipole moments in the 
quadrupole-removed maps. Gosmic variance dominates the uncertainty in the fitted 
a; consequently, we weight each pixel uniformly when computing the cross-correlation 
function. 

We test for possible biases in the cross-correlation technique by analyzing 2000 
additional simulations which include a signal T^''^ = 0.002 mK MJy"^ sr /diRbe 
addition to the GMB and noise. We recover a mean coefficient (a) = (1.96 ± 0.03) x 
10~^ mK MJy~^ sr at |6| > 20°, in excellent agreement with the simulation input. We 
have also repeated the analysis using the cross power spectra and obtained coefficients 
a in agreement with those found using the cross correlation function C{9). 
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3. Results 



Table 1 shows the fitted coefficient a from the cross correlation of the DMR 2-year 
maps with various Galactic template maps, compared to 2000 random realizations of 
CMB and instrument noise. We detect no statistically significant correlation between 
the DMR maps and any synchrotron template: the amplitude of any signal with the 
the 408 MHz survey or the cosmic-ray synchrotron model is smaller than the noise in 
the DMR maps. All of the DMR maps show a statistically significant correlation with 
the DIRBE far-infrared maps, indicating a detection of signal with common spatial 
structure in the high-latitude portion of the two data sets. 

The correlation coefficients are averages over the high-latitude sky and do 
not provide specific information on individual patches of the sky. Wc estimate the 
amplitude of the signal at each DMR frequency by multiplying the standard deviation 
of the Galactic template maps (after quadrupole subtraction) by the fitted coefficients 
a. Since the three DIRBE maps trace the same warm dust emission, the DIRBE 
templates provide nearly identical estimates of Galactic emission at each of the DMR 
frequencies. The uncertainties in the fitted coefficients a between any one DMR map 
and the three DIRBE maps are dominated by the noise and CMB signal in the DMR 
maps, and are not statistically independent. We therefore adopt the unweighted 
average of the fitted signal at each DMR frequency, with uncertainty estimated by 
the smallest uncertainty of the three DIRBE cross-correlations. 

Table 2 shows the estimated Galactic signal in the DMR maps traced by 
the two synchrotron templates and the DIRBE dust morphology. The Galactic 
signal correlated with the DIRBE dust distribution has rms fluctuations ATcai = 
7.1 ± 3.4 //K antenna temperature at |6| > 30° and 7° angular resolution in the most 
sensitive 53 GHz channels. For comparison, the rms instrument noise at 53 GHz is 
98 nK. Galactic emission is a minor contribution to the brightness flTictTiations at 
millimeter wavelengths. Specifying the template of the Galactic emission allows a 
significant increase in sensitivity to faint extended emission. 

Since the DIRBE maps are dominated by infrared dust emission, it is natural 
to interpret the correlated DMR-DIRBE signal as evidence for dust emission at 
millimeter wavelengths. However, the sharp rise in signal amplitude at 53 and 
31.5 GHz is inconsistent with emission by dust alone. The amplitude (in antenna 
temperature) of emission traced by the DIRBE maps more than doubles between 
53 and 31.5 GHz, requiring some emission component with negative spectral index. 
The spectral index T oc z/^ of the signal, without modelling any components, is 
/3 = -1.7 ± 0.4 between the 31.5 and 53 GHz data and (3 = -1.5 ± 0.9 between the 
53 and 90 GHz data (68% confidence), suggestive of an admixture of dust (/3 ^ +2) 
and free-free {(3 = —2.15) emission, but inconsistent with dust emission only. 

The DMR-DIRBE signals are not appreciably affected by synchrotron emission 
(/3synch ~ —3.1 at millimeter wavelengths). Direct cross-correlation of the DMR 
maps with either the 408 MHz survey or the cosmic-ray synchrotron model shows 
no statistically significant signal. The fitted coefficients for the 408 MHz survey at 
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|6| > 20° are significant only at 90% confidence, and have a spectrum closer to the 
CMB than to synchrotron emission (Table 2). The 408 MHz coefficients all become 
negative at the higher cut |6| > 30°, which we take as further evidence that the 
correlations at |6| > 20° are dominated by chance alignment of the CMB with a feature 
in the 408 MHz survey. Similarly, there is no compelling evidence for synchrotron 
emission in the 2-year DMR maps with spatial template traced by the cosmic-ray 
model. We adopt an upper limit for either model ATgynch < 12 //K (95% confidence) 
for fluctuations in synchrotron emission at 31.5 GHz. 

Cross-correlations between the DIRBE maps and the 408 MHz survey provide 
further evidence that the rising signal in the low-frequency DMR channels is not 
caused by synchrotron emission. We cross-correlate the DIRBE and 408 MHz 
template maps (with the 408 MHz survey scaled to 53 GHz using a spectral index 
P — —3 to provide direct comparison with the DMR 53 GHz results) and obtain 
a fitted coefficient at |6| > 30°, agynch = 3 x 10~^ mK MJy~^ sr, two orders of 
magnitude smaller than the correlation observed between the DMR and DIRBE maps. 
Synchrotron emission at high latitude does not strongly correlate with the DIRBE 
far-infrared dust emission. 

We conclude that synchrotron emission does not account for the observed 
correlation between the high-latitude DMR and DIRBE sky maps, leaving free-free 
emission as the most likely alternative. We decompose the sky signals from Table 
2 into dust and free-free components assuming /3 = 2 for the dust. Table 3 shows 
the rms dust and free-free normalization at 53 GHz. Approximately half of the 
high-latitude Galactic signal in the DMR 90 GHz data results from dust and half 
from free-free emission, while the 53 GHz data are dominated by free-free emission. 
Free-free emission increases markedly at the lower Galactic cut. Note that \b\ > 20° 
includes emission from the Orion and Ophiuchus complexes, which contain substantial 
amounts of ionized gas. The increase in fitted free-free amplitude when regions known 
to contain ionized gas are included in the analysis provides further support for the 
existence of free-free emission spatially correlated with the infrared dust emission. 

Tables 1 and 2 demonstrate a significant detection of dust and free- free emission 
in the high-latitude DMR data. We may use the combined DIRBE and DMR data 
to limit dust emission properties at long wavelengths. We fit the rms DIRBE and 
DMR signals to dust models of the form = t{^)^B^{T) +Ag{^)-^-^^ (model 1) or 
h = Ti{^fB,{Ti) + T2{^fB,,{T2) +Aff(i;^)-0-i^"(model2),i.e., a'model with a single 
dust population with enhanced submillimeter emissivity or a two-temperature model 
with emissivity. Both models include a free-free emission term Aq. Figure 1 shows 
the spectra of the intensity fiuctuations and the fitted dust and free-free emission for 
a single dust component (model 1) at |6| > 20°. Figure 2 shows the contours of this 
model in the temperature-emissivity plane. The best fit occurs for dust temperature 
T = 181? K and /3 = 1.91^;^ with r = (1.2l^j) x 10"^ (68% confidence), in agreement 
with the absolute FIRAS spectra (Wright et al. 1991, Reach et al. 1995a). The DMR 
data place a lower limit to a frequency-dependent dust emissivity over the wavelength 
range 6 mm to 100 /xm: /3>1.2 at 95% confidence. 
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The two-component dust model is not significantly constrained by the DMR- 
DIRBE cross-correlation. Provided one component is near 18 K, changes in the 
temperature of the cold component are compensated by correlated variations in 
opacity. We may fix the temperatures in model 2 to the values derived from the 
FIRAS |&|>30° spectra and compare the opacities derived from fluctuations in the 
warm and cold components to the total opacities derived from FIRAS (Table 4). The 
opacities for both the warm and cold components derived from the DMR-DIRBE 
cross-correlation are smaller than the corresponding values derived from the FIRAS 
spectra, as expected if the dust emission contains a significant monopole or quadrupole 
term. 

Figure 3 summarizes the contribution of various Galactic emission mechanisms 
(standard deviation of each component averaged over |6|>30°) at 7° angular 
resolution. The width of each band refiects the 68% confidence level uncertainty in 
the amplitude of each component, and includes the absolute calibration uncertainty 
of the template maps. The "cosmic-ray" synchrotron model provides a modestly 
better fit to the available low-frequency data than the 408 MHz survey scaled with 
a spatially invariant spectral index (Bennett et al. 1992). Accordingly, we normalize 
the synchrotron curve in Figure 3 by the cross-correlation of the DMR 31.5 GHz 
channel with the cosmic-ray model, as indicated by the filled circle. Fluctuations in 
the synchrotron emission at DMR frequencies are modestly smaller than expected; 
the cosmic-ray model predicts ATgynch = 13 /xK at |6|>30° in the 31.5 GHz map. We 
note that our results are insensitive to the zero level of the radio surveys, which can 
influence the amplitude of the predicted fluctuations through a monopole change in 
the spectral index. The 95% confidence level upper limit on synchrotron emission 
(ATsynch<12 /uK at 31.5 GHz) is the same for both the cosmic-ray and 408 MHz 
templates, and may be used to set an upper limit on the high-latitude synchrotron 
spectral index: /3synch< — 2.9 between 408 MHz and 31.5 GHz. 

Free-free emission is normalized by the free-free values in Table 2 and 
dominates fiuctuations in diffuse synchrotron emission for frequencies above 20 GHz. 
The dust spectrum inferred from the DMR-DIRBE cross-correlation is in agreement 
with more precise spectra from the COBE FIRAS experiment. The dust model in 
Figure 3 represents a range of high-latitude spectra from Tdust = 21 K and /3dust — 1-4 
to Tdust = 16 K and /9dust = 2, with opacity normalized to the spatial variation of the 
DIRBE 240 /im map. The upper limits on dust emission at millimeter wavelengths 
from the DMR-DIRBE cross-correlation dust model (after removing the free-free 
component) are also shown. 
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4. Discussion 



A significant result of this analysis is the detection of fluctuations in the microwave 
sky brightness with the angular variation of the DIRBE far-infrared maps and the 
frequency spectrum of free-free emission. The angular structure in the far-infrared 
maps is dominated by interstellar dust associated with atomic gas (Boulanger & 
Perault 1987). Our results suggest that the ionized component of the interstellar 
medium is at least partially correlated with the atomic gas. 

Previous comparisons of the Ha recombination line from the ionized gas 
revealed no clear correlation with the Hi hyperfine line from the neutral gas (Reynolds 
1987), although one filament detected in Ha may be related to a parallel offset 
structure in Hi (Reynolds et al. 1995). Similar correlated structures, as well as 
evaporating surfaces of clouds embedded in a hot ionized gas ( cf McKee & Ostriker 
1977) could account for the correlation we observe between the dust and free- 
free emission. The rms variation in the inferred free-free emission at 53 GHz is 
ATff = 5.3 ± 1.8 iiK at 7° angular resolution, averaged over the high-latitude sky 
(|6|>30°) after removal of the quadrupole component. The amplitude of this signal 
compares favorably with the predicted signal ATgf 7 /jK inferred from fluctuations 
in Ha emission (Reynolds 1992), where for convenience we reference ATg to 53 GHz. 

The detection of free-free emission correlated with the DIRBE far-infrared 
maps provides a template for at least one component of Galactic free-free emission. 
We may place an upper limit for the combined free-free emission from all sources 
by forming a linear combination of the DMR maps designed to pass emission with 
spectral index -2.15, cancel emission with a CMB spectrum, and minimize instrument 
noise: 

Tff = 0.40 X i(r^,A ± T^^b) - 0.09 X ^(T^^j, ± T^,^) - 0.37 x ^(T^,^ ± T^,^), 

where T' is the antenna temperature in each DMR channel after subtracting 
synchrotron and dust emission using the cosmic-ray and DIRBE models, respectively. 
Since both synchrotron and dust are small compared to the free-free emission at 53 
GHz (Figure 3), their residual contribution after correction is negligible. This linear 
combination has only 33% larger noise than the most sensitive 53 GHz channels. We 
smooth the maps with a 7° FWHM Gaussian to reduce further the effects of noise, 
remove a fitted quadrupole, and compare the variance of the (A+B)/2 sum map 
to the (A-B)/2 difference map to obtain an estimate for the free-free fluctuations 
ATff = 3.9 ± 5.7 nK, in excellent agreement with the value ATg = 5.1 ± 1.6 /xK from 
the DMR-DIRBE cross-correlation at the same 10° effective smoothing. Free-free 
emission from sources uncorrelated with the far-infrared dust emission, although 
present at high latitudes, does not dominate the component correlated with the dust 
emission for the range of angular scales probed here. 

Figure 4a shows the power spectrum of the inferred free-free emission at 53 
GHz. We calculate the power spectrum of the high-latitude DIRBE 240 //m map 
using a set of basis functions orthonormal on the region |6|>30° (Gorski 1994), and 
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normalize to free-free emission at 53 GHz using the correlation coefficient derived 
from a dust /free- free decomposition of the correlation coefficients from Table 1, 
ciflf = (3.35 ± 1.11) X 10~^ mK MJy~^ sr. The power falls rapidly at small angular 
scales: 



where P{£) — 2;^J2m(4m mean power at multipole order i {£ n/O), and 

the quadrupole normalization 02 = (6.02 ± 3.99) /xK^. The power spectrum in Figure 
4 is an average over the high-latitude sky. On smaller angular scales, Gautier et al. 
(1992) estimate the dust power spectrum using 8° x 8° IRAS patches and find P oc 
as well. 

The DMR and DIRBE data provide the normalization for the correlated free- 
free emission on angular scales larger than 7°. The dust power spectrum continues 
to smaller angular scales, and it is plausible that the spatial correlation between the 
dust and the warm ionized gas also persists. We thus estimate the rms free-free 
fluctuations on various angular scales by integrating the free-free power spectrum 
in Figure 4a over the range in i corresponding to the patch size and the Gaussian 
beam width. Table 5 compares the model predictions to Ha and radio measurements. 
An Ha map of a single 12° x 10° patch suggests rms variations ATg = 1.4 /iK at 
0?8 angular resolution (Reynolds 1992). Similar observations of a 7° x 7° patch at 
0?1 angular resolution yield an upper limit of ATfj<1.2 /xK (Gaustad et al. 1995), 
while deep VLA observations at 1' resolution provide an upper limit ATff<1.3 fiK 
(Fomalont et al. 1993). There is generally good agreement between observations and 
the power spectrum predictions, suggesting that the dominant free-free emission is 
correlated with the dust and that the correlation persists to smaller angular scales. 

If the free-free/dust correlation docs persist to smaller angular scales, it 
has interesting impfications for experiments measuring CMB anisotropy at degree 
angular scales. Measurements of CMB anisotropy must contend with the problem 
of estimating the foreground Galactic emission and removing it if it contributes 
a significant fraction of the CMB anisotropy. Since no template for the free-free 
emission exists, experimenters typically observe in several frequency bands and use a 
linear combination of the angular anisotropy at different frequencies to remove any 
emission with the free-free emission spectrum. Although this removes the foreground 
emission, it increases the noise in the corrected CMB results, often to the point 
of obscuring the CMB signal. If free-free emission could be demonstrated to be a 
negligible contaminant, significant reductions in noise could be realized. 

Table 5 shows the free-free power spectrum predictions for two cases of interest. 
A full-sky map at 53 GHz with angular resolution 0?5 would observe rms free- 
free emission of 6.1 ± 2.0 fiK averaged over the high- latitude sky after quadrupole 
subtraction. Surveys restricted to smaller regions exclude power from Fourier modes 
with wavelength larger than the survey size. A 10° x 10° survey (corresponding 
to modes £>18) with angular resolution 0?5 would measure rms free-free emission 
of 1.5 ± 0.5 //K, where the quoted errors reflect only the uncertainty in the free- 
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free normalization and do not include the variation from patch to patch on the sky. 
Emission at this level is well below the expected magnitude of CMB anisotropy, 
although individual patches may show free-free emission significantly larger than the 
high-latitude average. 

The small amplitude of the free-free signal explains why no structure is 
observed in the maps of free-free emission derived from linear combinations of the 
two-year DMR maps. The instrument noise in the free-free map presented above 
is 20 times larger than the free-free signal estimated from the DMR-DIRBE cross- 
correlation. Figure 4b shows the power spectrum of this free-free map, including the 
uncertainty from instrument noise. It is clear that the noise per multipole order I is 
much larger than the signal. The increase in sensitivity required to detect signals at 
the few yuK level results from the a priori specification of the DIRBE map as a spatial 
template for the emission. 

The high-latitude Galactic signal is small compared to the CMB anisotropy 
in the DMR maps and does not significantly alter previously published estimates 
based on the uncorrected DMR 53 GHz sky maps. The standard deviation at 
7° angular resolution of the combined GMB and Galactic signals at |6|>30° is ATgky = 
42.6 ± 8.4 //K antenna temperature at 53 GHz. The CMB and Galactic signals are 
uncorrelated; correcting ATgky for the estimated Galactic signal ATcai = 7.1 ±3.4 /iK 
yields an estimated CMB anisotropy ATcmb = 42.0 ± 8.6 ;uK, within 0.7 /xK (7% 
of the statistical uncertainty) of the uncorrected result. The CMB anisotropy may 
also be characterized by a power-law distribution of primordial density perturbations 
P{k) oc /c" with quadrupole normalization Qrms-ps ^-nd index n. Removing the 
Galactic emission traced by the DIRBE maps changes Qrms-PS or n by less than 0.3 
standard deviations, well within their respective statistical uncertainties. 



5. Conclusions 

We detect statistically significant cross-correlation between the COBE iwoycai DMR 
maps and the DIRBE 240, 140, and 100 ^m. maps. The rms amplitude of the 
correlated Galactic signal in the DMR maps at |6| > 30°, after quadrupole subtraction, 
is 15.6±6.0, 7.1±3.4, and 3.9±3.8 jiK at 31.5, 53, and 90 GHz, corresponding to dust 
emission ATdust = 0.9 ± 1.3 /iK and free-free emission ATff = 5.3 ± 1.8 jiK at 53 GHz 
and 7° angular resolution. The dust emission may be used to put a lower limit on the 
emissivity at millimeter wavelengths, AT^ust oc with {3 > 1.2 at 95% confidence. 
The dust temperature and emissivity derived from fluctuations in the DMR and 
DIRBE maps are in agreement with values derived from the absolute FIRAS spectra 
at wavelengths between the DMR and DIRBE coverage. Similar analyses using the 
408 MHz survey or a model of synchrotron emission based on the 408 MHz survey 
with spatially varying, frequency-dependent spectral index yielded upper limits for 
fluctuations in synchrotron emission of ATgynch < 12 //K at 31.5 GHz, corresponding 
to a spectral index /3synch < —2.9 between 408 MHz and 31.5 GHz. 
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We detect a component of Galactic free-free emission correlated with the 
DIRBE maps, providing a template for the angular distribution of free-free emission 
on angular scales larger than the 7° DMR beam. The rms variations in this correlated 
component, ATg = 5.3 ± 1.8 /xK, are in good agreement with the rms variations from 
all sources of free-free emission, ATgf = 3.9±5.7 //K, derived from a linear combination 
of the DMR 2-year maps. Extrapolating the power spectrum of the free- free emission 
to smaller angular scales yields an estimated contribution for rms fluctuations at 
0?5 angular scale observed at 53 GHz of AT^ = 6.1 ± 2.0 /xK for a full-sky map at 
\b\ > 30°, or ATff = 1.5 ± 0.5 /iK for observations restricted to a 10° x 10° patch. 

If the spatial correlation between the dust and warm ionized gas observed on 
large angular scales persists to smaller angular scales, free-free emission should not 
be a serious contaminant to measurements of the medium-scale CMB anisotropy. It 
is worth pointing out, however, that small regions of the sky have variations much 
larger than the standard deviation over the full sky. If free-free emission is correlated 
with dust emission on all angular scales of interest, the IRAS or DIRBE maps may 
be used to estimate both dust and free-free emission in individual regions of the sky. 
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Table 1: DMR-Galactic Template Cross- Correlation Coefficients" 



Template Map DMR Map 10^ a (|&|>20°) 



408 MHz" 



10^ a (|6|>30°) 



31.5 GHz 

53 GHz 
90 GHz 



2.32 ± 1.57 
1.54 ± 0.91 
1.14 ± 1.00 



-0.70 ± 2.25 

-0.28 ± 1.31 
-0.068 ± 1.43 



Cosmic-Ray'^ 



DIRBE 240 nm!^ 



DIRBE 140 ^xm.'^ 



DIRBE 100 //m'* 



31.5 


GHz 


(0.21 ± 0.40) X 


10^ 


(0.63 ± 0.46) X 


10^ 


53 


GHz 


(0.13 ± 0.93) X 


10^ 


(0.78 ± 1.07) X 


10^ 


90 


GHz 


(3.83 ± 5.25) X 


10^ 


(-2.24 ± 5.95) X 


10^ 


31.5 


GHz 


12.63 ± 1.80 




10.09 ± 3.65 




53 


GHz 


5.61 ± 1.03 




4.03 ± 2.06 




90 


GHz 


2.34 ± 1.20 




2.34 ± 2.28 




31.5 


GHz 


9.27 ± 1.33 




7.38 ± 2.90 




53 


GHz 


4.06 ± 0.77 




3.90 ± 1.70 




90 


GHz 


1.82 ± 0.86 




1.47 ± 1.83 




31.5 


GHz 


18.06 ± 2.54 




14.57 ± 6.03 




53 


GHz 


6.88 ± 1.40 




6.46 ± 3.45 




90 


GHz 


2.76 ± 1.61 




4.56 ± 3.89 





" The antenna temperature at each DMR frequency of Galactic emission traced by 

template map X is T4 = aX mK (see text). 

^ a has units mK since the template map has units K. 

a is dimensionless since the template map has units mK. 
^ a has units mK (MJy/sr)~^ since the template map has units MJy sr~^. 
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Tabic 2: RMS Galactic Signal in DMR Sky Maps 



Template Map 


DMR Frequency 
(GHz) 


rms at |6|>20° 
(/xK) 


rms at |6|>30° 
(/iK) 


408 MHz 


31.5 


7.6 ± 5.1 


-1.9 ± 6.1 




53 


5.0 ± 3.0 


-0.8 ± 3.6 




90 


3.7 ± 3.3 


-1.9 ± 3.9 


Cosmic-Ray 


31.5 


2.9 ± 5.4 


8.0 ± 5.8 




53 


0.3 ± 2.8 


2.2 ± 3.0 




90 


2.4 ± 3.2 


-1.3 ± 3.5 


DIRBE 


31.5 


35.2 ± 4.9 


15.6 ± 6.0 




53 


14.8 ± 2.7 


7.1 ± 3.4 




90 


6.3 ± 3.1 


3.9 ± 3.8 



Table 3: Pree-Pree and Dust Antenna Temperature at 53 GHz 



Component 


rms at |6|>20° 


rms at |6|>30° 






(/.K) 


Pree-Pree 


12.0 ± 1.5 


5.3 ± 1.8 


Dust 


1.1 ± 1.1 


0.9 ± 1.3 



Table 4: Dust Opacities in Two-Component Model, |6|>30' 



Data Set 


Warm Opacity 


Cold Opacity 




10^ Ti 


10^ T2 


FIRAS" 


2.18 ± 0.03 


15.0 ± 0.7 


DMR-DIRBE^ 


0.87 ± 0.08 


2-3^2.3 



" PIRAS values derived from mean dust spectra (Reach et al. 1995a). 
DMR-DIRBE values derived from intensity fluctuations. 
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Table 5: Free-Free Observations and Model Predictions' 



Angular Resolution 


Patch Size 


Observed ATfj 


Predicted AT^ 


Reference 


10° 


|6|>30° 


3.9 ± 5.7 


5.3 ± 1.6 


This work 


0?8 


12° X 10° 


1.4 


1.5 ± 0.5 


Reynolds 1992 


0?1 


7° X 7° 


< 1.2 


1.4 ± 0.5 


Gaustad et al. 1995 


1' 


7' X 7' 


< 1.3 


0.15 ± 0.05 


Fomalont et al. 1993 


0?5 


10° X 10° 




1.5 ± 0.5 




0?5 


|6|>30° 




6.1 ± 2.0 





" All free-free values are in antenna temperature at 53 GHz 

rms fluctuations estimated from power spectrum of free-free emission correlated 
with DIRBE dust emission (see text). 
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i30 3 0.3 0.06 




10 100 1000 5000 

Frequency (GHz) 

Figure 1: Spectrum of correlated intensity fluctuations in the DMR and DIRBE 
data at |6|>20°. Far-infrared points represent the standard deviation of the DIRBE 
maps after subtracting a fitted monopole, dipole, and quadrupole, and include 
calibration uncertainties. Long-wavelcngth points are the values inferred from the 
cross-correlation of the DMR and DIRBE maps (Table 2). The fitted free-free 
emission and single- component dust model are also shown. 
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Figure 2: contours for the temperature and emissivity index of a single dust 
component fitted to the DMR-DIRBE intensity fluctuations at |6|>20°. The asterisk 
marks the best fitted value. Free-free emission is fitted separately for each value of 
the dust emissivity. 
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Figure 3: Frequency spectra of high-latitude Galactic emission {rms fluctuations at 
|6|>30° after quadrupole subtraction for a 7° FWHM beam). Free-free emission 
dominates synchrotron emission at frequencies above 20 GHz. Upper limits on dust 
emission at millimeter wavelengths from the DMR-DIRBE cross-correlation, after 
correction for free-free emission, are also shown. 
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Figure 4: Mean power spectrum of the inferred free-free emission at 53 GHz. (top) 
Free-free power spectrum derived from the DMR-DIRBE cross-correlation at Galactic 
latitude |6|>30°. The filled circles show the DIRBE power spectrum on angular scales 
larger than 7°, normalized to the correlation coefficient of the free-free emission at 53 
GHz (see text). The sohd hne shows a fitted spectrum, (bottom) Mean power 
spectrum of a free-free map derived from a linear combination of DMR maps without 
specifying the DIRBE maps as the spatial template. Open circles are the actual 
power at each multipole order, while the grey band shows the 68% confidence interval 
from instrument noise. Filled circles are26he power spectrum of the DMR-DIRBE 



